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PREFACE 


The  work  described  in  this  report  on  studies  of  a  laser/nuclear  thermal  hardened 
body  armor  was  undertaken  during  the  period  January  31,  1991  to  September  30,  1991. 

The  professional  affiliation  of  N.Y.  Misconi  is  Center  for  Geo-Space  Environmental 
Research,  Department  of  Mechanical  &  Aerospace  Engineering,  Florida  Institute  of 
Technology,  Melbourne,  Florida.  The  Natick  affiliation  of  Gerald  J.  Caldarella  and 
Joseph  F.  Roach  is  Physics  and  Engineering  Branch,  Fiber  ai.d  Polymer  Division,  Soldier 
Science  Directorate. 

The  funding  for  this  research  was  Program  Element  62786,  Project  Number  ILIAH98, 
Task  Number  CAB0O. 

The  citation  of  tradenames  in  this  report  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  an  item. 

The  authors  are  grateful  to  Marcia  Lightbody  for  her  assistance  in  editing  and 
preparing  this  document. 


STUDIES  OF  A  LASER/NUCLEAR  THERMAL-HARDENED  BODY  ARMOR 


SUMMARY 

This  final  report  contains  three  parts.  Part  1  gives  the  history  and  background 
concerning  the  earlier  research  supported  by  the  Air  Force  Office  of  Scientific  Research 
(AFOSR),  during  the  years  1984  to  1989.  This  research  was  under  AFOSR’s  Satellite 
Survivability  Program,  which  aimed  at  developing  a  laser  shield  to  harden  space  systems 
against  space-  and  ground-based  laser  weapons.  We  used  the  concept  of  particle  light 
scattering  instead  of  particle  light  absorption  in  developing  this  shield.  The  light 
scattering  concept  utilized  the  makeup  of  a  layer  of  highly  irregular  pm  sized 
particles  that  are  highly  pure.  The  multiple  scattering  insVe  the  layer  by  the  photons 
enables  the  majority  of  them  (-99%)  to  be  reflected  back  and  away  from  the  target.  An  18 
watt,  continuous  w'ave  (CW)  argon  laser  was  used  in  this  investigation.  This  laser  shield 
protects  targets  against  laser  weapons  from  the  0.22  to  2.4pm  region  of  the 
spectrum. 

Part  11  of  the  report  deals  with  protecting  targets  against  the  CO-,  laser 
(10.6  pm  )  by  utilizing  the  same  concepts  mentioned  above.  The  layers  were 
made  out  of  naturally  Recurring  NaCl  particles,  which  have  very  low  absorption 
coefficient  (-10'  cm'  )  in  the  region  of  the  spectrum  between  2  to 
16  pm  .  The  natural  salt  particle  layers  were  subjected  to  a  COi,  20 
watt,  CW  laser  and  found  to  reject  more  than  -99%  of  the  incoming  radiation. 

Part  111  of  '.his  report  contains  the  details  of  a  created  mathematical  model  based 
on  radiative  transfer  equations  to  calculate  the  temperature  profile  of  the  layers  as 
they  are  subjected  to  incoming  radiation.  A  computer  code  was  developed  to  perform  these 
model  calculations,  in  FORTRAN,  which  is  included  in  She  report.  The  computer  code 
accepts  input  parameters,  such  as  particle  sizes,  refractive  index,  absorption 
coefficient,  etc.  It  outputs  parameters,  such  as  the  amount  of  radiation  that  emerges 
out  of  the  layers,  temperature  distribution  across  the  layer,  melting  threshold,  etc. 

The  purpose  of  the  model  is  to  predict  the  best  combinations  of  parameters  to  optimize 
the  radiation  by  these  layers. 
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Part  I.  RESEARCH  ON  THE  INTERACTION  OF  SOLID  PARTICLES 
WITH  USER  BEAMS 


INTRODUCTION  AND  HISTORY 

From  1985  to  1989  funds  from  the  research  program  "The  Interaction  of  Solid 
Particles  with  Laser  Beams:  Application  to  the  Defense  of  Satellites",  under  "The  Satellite 
Survivability  Program"  (SSP)  at  the  Air  Force  Office  of  Scientific  Research  (AFOSR)  were 
used  to  build  a  state-of-the-art  Laser -Par tide  Dynamics  facility  at  the  University  of 
Florida.  The  research  resulted  in  a  patent  that  is  pending  for  the  material,  process,  and 
later  products  that  involve  protecting  many  types  of  military  targets  against  lethal  laser 
weapons. 

Since  the  onset  of  research  for  finding  methods  to  protect  satellites  against  laser 
weapons,  absorption  of  the  laser  energy  before  reaching  the  target  has  been  the  prevailing 
approach.  Our  research,  however,  was  to  find  a  way  tc  scatter  the  laser  beam,  so  that,  at 
most,  only  a  very  small  percentage  could  reach  the  satellite,  thus  rendering  the  weapon 
harmless.  This  approach  has  the  advantage  that  the  shield  is  not  damaged  in  the  process  of 
protecting  the  target,  and  thus  is  reusable.  Although  the  main  thrust  of  our  research 
would  contribute  to  the  understanding  of  scattering  by  a  cloud  of  particles,  we  also 
studied  closed  aggregates  of  particles,  in  the  form  of  a  "layer",  many  particles  thick  and 
flattened  on  both  top  and  bottom  surfaces.  Our  research  on  layer  scattering,  to  date,  has 
emphasized  the  distribution  of  scattered  light  as  a  function  of  particle  size,  layer 
thickness,  and  material.  It  is  this  layer  scattering  facet  of  our  research  that  was  shown 
to  be  the  most  promising  in  the  protection  of  targets  from  laser  beams. 

A  thin  layer  of  small,  highly  transparent  particles  will  act  like  a  diffuse  mirror 
when  exposed  to  an  incident  laser  beam  (see  Annual  Report  1  -  AFOSR  Contract 
F49620-85-C-01 17).  The  low  absorptivity  of  the  qiaterial  prevents  damage  to  the  layer,  even 
when  exposed  to  intensities  as  high  as  1.5  MW/cm4'.  We  have  found  it  straightforward  to 
create  layers  no  more  than  1  mm  thick  that  will  scatter  more  than  99%  of  the  laser 
radiation  back  in  the  direction  of  the  incident  beam.  The  thinness  of  the  layers  and  their 
porosity  make  for  a  very  lightweight  shield  that  would  be  desirable  for  human  body  armor. 

The  following  sections  describe  the  nature  of  light  scattering  by  a  single  particle 
and  by  many  particles  packed  in  a  layer.  In  the  former,  we  try  to  reproduce  the  Mie 
scattering  using  laser-particle  levitation  technique.  In  the  latter  we  measure  the  angular 
distribution  of  the  scattered  light  above  and  below  a  well-packed  layer  of  highly  irregular 
silica  particles.  These  measurements  are  then  used  to  calculate  the  laser  beam  suppression 
ratio  by  the  layer.  This  ratio  is  the  light  intensity  measured  below  the  layer 
(backscattering  of  the  laser  beam)  divided  by  the  light  intensity  above  the  layer. 

Comparisons  of  this  ratio  are  made  between  several  types  of  silica  particles,  different 
particle  size  range,  and  different  laver  thickness.  The  purpose  of  these  comparisons  is  to 
optimize  the  suppression  of  the  laser  beam, 


THEORETICAL  APPROACH  TO  L1CIIT  SCATTERING 

Theorists  still  have  not  been  able  to  solve  the  problem  of  light  scattering  by  even 
a  single  irregularly  shaped  particle.  In  fact,  exact  solutions  are  only  available  for 
single  spheres  and  small  spheroids,  and  even  these  require  time-consuming  computations.  To 
accommodate  effects  such  as  interference  and  diffraction,  one  must  sum  the  complex 
amplitudes  of  each  ray  entering  the  layer,  taking  into  account  the  phase  shift  across 
adjacent  paths,  not  just  sum  the  intensity  of  the  beam.  With  the  huge  number  of  randomly 


oriented  irregular  particles,  even  today’s  best  supercomputers  are  not  sufficient  to 
give  an  exact  solution  for  multiple  scattering  in  a  layer. 

Figure  1  shows  the  theoretical  Mie  Scattering  curve  of  a  highly  pure  silica 
sphere  of  33.0  jum  and  refractive  index  n  =  1.496,  superimposed  on  the 
experimental  results  found  from  our  laser-levitation  experiment  .  Laser-levitation 
is  a  technique  used  at  FIT  to  study  light  scattering  from  single  particles  and  sometimes 
doublets  or  triplets.  The  agreement  between  the  experimental  values  and  the  theoretical 
curve  is  truly  remarkable. 

The  problem  of  making  theoretical  scattering  functions  for  an  ensemble  of 
particles  or  a  layer  of  particles  can  not  be  solved  even  with  today’s  fastest 
supercomputers.  To  demonstrate  the  difficulty  of  this  problem,  consider  the  scattering 
of  sunlight  from  interplanetary  dust  particles  (otherwise  known  as  the  zodiacal  light). 
Here  we  have  scattering  of  sunlight  from  a  cloud  of  irregularly  shaped  particles,  and 
the  only  way  to  extract  information  on  scattering  by  these  particles  is  by  inverting  the 
brightness  integral  to  obtain  an  empirical  scattering  function  for  the  cloud.  Figure  2 
shows  the  different  empirical  scattering  curves  obtained  by  different  observers  using 
the  previously  mentioned  method.  It  is  obvious  that  many  features  of  the  Mie  curve 
(Fig.l)  are  totally  washed  out  in  Figure  2.  There  is  even  doubt  on  the  validity  of 
obtaining  an  emperical  scattering  function,  since  by  inverting  the  brightness  integral, 
one  has  to  make  assumptions  on  the  particles  density  distribution  and  size.  Considering 
the  above,  we  are  resigned  to  the  fact  that  the  distribution  of  light  scattered  from  a 
layer  of  particles  can  be  found  only  by  measuring  it  experimentally. 


Figure  1.  Scattering  of  514.5  nm  laser  light  from  a  33  jura  glass 
sphere.  Measurements  (dots)  are  compared  to  theory  (solid  line). 
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Figure  2.  Zodiacal  light  scattering  curves.2*^ 

RESULTS:  INTERACTION  OF  SILICATE  PARTICLES  WITH  LASER  BEAMS 

We  tested  many  different  layers  of  small  particles.  Each  layer  was  formed  from 
some  variety  of  silicate,  be  it  glass  or  silicate  sand.  The  sand  samples  were  rinsed  in 
HNOj  and  then  distilled  water  to  remove  any  brine  residues.  The  material  was  crushed 
to  a  nne  powder,  then  sorted  into  various  particle  size  ranges  using  sieves  from  250 
down  to  6S/m  The  particles  were  then  packed  into  wells  built  on  microscope 
slides,  w-ith  depths  of  0.25,  0,50.  I.  and  2  millimeters,  to  give  layers  of  these 
thicknesses. 

Particles  less  than  lOMm  frequently  clung  to  the  larger,  sieved 
particles.  W'e  "cleaned"  some  of  these  samples  by  rinsing  with  distilled  water,  and  left 
others  "dirty"  to  test  the  difference. 

A  silicon  photodiode  detector  w-as  mounted  on  a  goniometer  arm  and  centered  in 
the  layer  as  in  Figure  3.  The  beam  of  a  20  W  argon  ion  laser  was  brought  from  the 
bottom  side  of  the  sample,  as  shown  in  the  figure,  and  passed  through  the  glass  slide 
before  entering  the  layer.  W‘e  retested  several  samples  with  the  beam  incident  upon  the 
sample  from  the  top.  so  that  the  beam  was  intercepted  by  the  layer  before  reaching  the 
slide.  There  was  no  difference  in  the  readings. 
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Figure  3.  Experimental  setup  used  to  measure  layer  scattering 

Readings  were  taken  by  moving  the  goniometer  through  scattering  angles  (6) 
from  -165°  to  *165°  (6-0  indicated  transmission  in  the  direction  of  the  beam). 

The  readings  were  limited  at  the  high  scattering  angles  by  the  blocking  of  the  beam  by 
the  detector.  The  scattering,  as  expected,  was  symmetric  in  *6  .  However,  when 
the  layer  was  tilted  we  found  that  the  scattering  curve  became  symmetric  about  a  line 
perpendicular  to  the  layer,  and  not  about  the  beam  axis.  An  example  of  a  scattering 
curve  is  shown  in  Figure  4a.  When  corrected  for  the  projection  of  the  surface  as  seen 
by  the  detector,  the  angular  distribution  of  Figure  4b  is  obtained.  Although  the 
magnitudes  of  the  transmitted  and  reflected  hemispheres  varied  with  particle  size  and 
material,  the  general  profile  of  all  of  the  curves  was  functionally  similar  to  these  in 
Figure  4. 

The  scattering  curve  was  extrapolated  in  the  region  between  165°  and  180°  and 
then  integrated  to  determine  the  total  amount  of  light  reflected  vs.  the  total  amount 
transmitted.  This  was  then  reduced  to  a  beam  suppression  ratio  (BSR),  where  higher 
numbers  refer  to  a  larger  amount  of  light  reflected,  and  thus  better  protection  of  a 
target.  Values  above  100  were  found,  but  as  the  BSR  increased,  the  error  in  measurement 
also  increased,  due  to  the  sensitivity  limit  of  our  detector.  It  is  conceivable  that 
BSR  values  greater  than  10,000  were  achieved,  although  transmission  of  such  small  light 
levels  could  not  be  measured  with  the  apparatus  used  then. 
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Figure  4a.  Scattering  and  transmission  of  a  1  mm  thick  layer  of  Suprasil" 
particles  in  the  size  range  90-]25jum  . 


Beam  Suppression  for  Different  Types  of  Glass 

Figure  5  shows  our  experimental  results  on  beam  suppression  ratios  of  different 
kinds  of  glass  and  for  different  layer  thicknesses.  It  leaves  no  doubt  that  natural 
beach  sand  can  give  BSR’s  much  better  than  any  man-made  glass.  It  is  our  contention 
here  that  the  beach  sand  gives  a  higher  value  of  the  BSR  because  it  has  crystalline 
structure  that  is  not  destroyed  in  the  manufacturing  process,  not  because  of  an  increase 
in  absorption.  This  structure  results  in  a  better  light  reflection  capability  for  the 
sand  over  Suprasil™  and  the  others.  This  contention  is  supported  by  an  X-ray 
spectroscopy  analysis  performed  at  the  Department  of  Materials  Science  and  Engineering 
of  the  University  of  Florida,  where  measurements  of  various  samples  show  that  some  sand 
samples  are  more  pure  than  even  Suprasil  glass. 

To  test  the  level  of  protection  afforded  by  these  materials,  and  indirectly, 
their  absorption,  we  focused  the  full  power  of  the  laser  beam  just  below  the  surface  of 
each  sample,  giving  a  power  density  of  1.5  MW/cm2  of  continuous  energy  for  three 
minutes.  Most  samples  melted  under  this  power,  but  some  survived,  including  Suprasil 
glass,  which  we  knew  to  have  a  low  absorption  coefficient.  We  compared  the  scattering 
curves  of  Suprasil  with  others  that  also  were  undamaged,  and  found  that  the  naturally 
occurring  sands,  which  reflected  more  light  than  Suprasil  to  achieve  higher  BSR’s, 
corresponded  to  these  undamaged  samples. 


Naiuraf  Sana  Glass  Spheres 

Figure  5.  Beam  suppression  ratios  for  various  materials. 
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Particle  Shape 


Our  experimental  results  for  the  effect  of  particle  shape  on  the  fraction  of  the 
laser  beam  that  emerges  from  the  layer  in  the  forward  direction  are  shown  in  Figure  6 
Layers  of  the  same  thickness  but  made  up  of  irregular  particles  give  BSRs  a  factor  of 
seven  higher  than  equivalently  sized  spheres.  The  irregulars  reflect  the  laser  beam 
significantly  better  than  spheres. 


Figure  6.  Effect  of  particle  shape  and  laser  beam  suppression. 


Wavelength  Dependence 

Our  experience  with  glasses  and  silicates  showed  no  wavelength  dependence  (to 
v  ithin  the  accuracy  of  our  measurements)  across  the  .458  to  .514  4m  range  of 
our  test  laser.  This,  however,  is  a  small  fraction  of  the  .22  to  2.2  4m 
range  of  transparency  for  silicates.  Theoretically,  this  material  should  show  a 
functional  dependence  on  wavelength  due  to  a  change  in  the  particle  size  to  wavelength 
ratio,  as  well  as  because  of  a  wavelength  dependence  in  the  value  of  the  index  of 
refraction.  This  conjecture  needs  to  be  demonstrated  experimentally. 


Particle  Size  Dependence 

Our  preliminary  work  demonstrated  a  strong  inverse  relationship  between  the  BSR 
and  particle  size,  which  increased  more  rapidly  and  nonlinearly  toward  the  small  end  of 
our  particle  size  spectrum.  The  smallest  size  sieve  used  in  our  sizing  process  was 
684m,  This  meant  that  our  smallest  size  range  included  particles  from  less 
than  14m  up  to  684m,  While  we  know-  that  the  optimum  size  range  is 
less  than  684m,  more  experiments  need  to  be  performed  to  determine  if  there 
is  a  lower  limit. 
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Layer  Thickness 


More  than  any  other  parameter,  much  is  already  known  about  the  relationship  of 
layer  thickness  and  beam  suppression.  In  Figure  7  we  show  the  beam  suppression  ratios 
of  Suprasil  as  a  function  of  layer  thicknesses  and  also  as  a  function  of  different 
sizes.  Similarly,  in  Figure  8  we  show  the  same  dependence  except  for  natural  sand. 
Obviously,  as  thickness  increases,  less  light  is  transmitted.  The  one  big  problem  which 
remains  is  to  relate  thickness  to  absolute  numbers  for  the  cases  studied  above.  For 
example,  just  how  thick  does  a  layer  need  to  be  to  give  a  BSR  of  10,  or  of  100,  for  any 
given  material  and  size  range. 

One  other  question  is  of  interest:  does  layer  thickness  affect  melting,  and  if  • 
so,  how?  The  total  absorption  of  a  layer  depends  on  the  thickness  of  the  layer. 

However,  as  thickness  increases,  the  total  energy  absorbed  is  distributed  among  more  and 
more  particles,  and  the  trade-off  between  these  two  effects  is  not  clear.  Measurements 
are  needed  to  determine  the  effect  of  layer  thickness  on  melting. 


PARTICLE  SIZE  (pm) 


Figure  7.  BSR  as  a  function  of  particle  size  for  layers  of  Suprasil  glass. 
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Figure  8.  BSR  as  a  function  of  particle  size  for  natural  sand. 


Part  II.  RESEARCH  ON  LASER/NUCLEA^  THERMAL-HARDENED 

BODY  ARMOR 


INTRODUCTION 

This  research  project  is  designed  to  extend  our  prior  research  on  scattering  by  a 
layer  of  small  (<  250ym  )  particles  to  include  NaCl  particles.  A  crushed 
layer  of  NaCl  particles  is  subjected  to  the  radiation  from  a  20  watt  continuous  wave 
(CW)  C02  laser  to  determine  the  suppression  of  this  radiation  for  the  purpose  of  using 
these  particles  in  developing  a  laser/nuclear  thermal  hardened  body  armor.  Tests  are 
also  made  to  determine  the  melting  thresholds  of  these  particles  under  intense 
radiation.  A  mathematical  model  appears  in  Part  111  that  was  developed  to  optimize  the 
reflectivity  by  using  these  particles  in  various  size  ranges  near  the  wavelength  of  the 
incident  laser  beam.  Prior  experiments  using  highly  pu;e  silica  particles  indicate  that 
we  can  reject  laser  light  to  such  a  degree  that  less  than  one  part  in  10  reaches  the 
target. 


The  study  involves  making  measurements  across  a  wide  dynamic  range.  Namely,  to 
accurately  access  the  effects  of  particle  size,  layer  thickness,  index  of  refraction, 
and  wavelength  on  the  distribution  of  transmitted  and  reflected  light  and  absorption  for 
very  low  transmission  levels.  Our  ultimate  goal  is  to  create  a  suitable  method 
Utilizing  this  unique  multiple  scattering  concept  to  improve  the  present  hardened  body 
armor  against  nuclear  and  laser  threats. 


RESEARCH  METHODOLOGY 

A  detector  for  very  low  energy  levels,  10*^  W,  was  mounted  on  a  goniometer 
arm,  centered  in  the  layer  as  in  the  sketch  of  Figure  3.  The  beam  of  a  20  watt  CW  C02 
laser  was  brought  from  the  bottom  side  of  the  sample,  as  shown  in  Figure  3,  and  passed 
through  the  zinc  selenium  (ZnSe)  slide  before  entering  the  layer. 

The  C02  laser  was  pulsed  to  avoid  damage  to  the  detector  (Pyroelectric' 
with  chopping  capability).  The  parameters  used  in  this  computer  pulsing  are  given  in 
Table  I. 


TABLE  1 

Power  Computations  of  C02  Laser.  Constant  Pulse  Repetition  Interval  is  100  ms 


Duty  Cycle 

Power 

% 

mW 

i 

200 

2 

400 

3 

600 

(Continued) 
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Table  1  (Continued) 


4 

800 

5 

1000 

6 

1200 

7 

1400 

8 

1600 

9 

1800 

10 

2000 

12 

2400 

14 

2800 

1  6 

3200 

18 

3600 

20 

4000 

22 

4400 

24 

4800 

26 

5200 

28 

5600 

38 

7600 

48 

9600 

58 

11600 

68 

13600 

78 

15600 

88 

17600 

98 

19600 

100 

20000 

RESULTS:  INTERACTION  OF  NiCI  PARTICLES  WITH  C02  LASER 

We  used  several  thicknesses  of  the  layer  of  particles  size  25  Mm 
starting  with  thickness  of  4  mm.  We  could  not  detect  any  radiation  coming  out  of  the 
layer  and  into  the  detector  Le^the  amount  of  radiation  reaching  the  detector  is  in  the 
noise  of  the  detector  (<  10'"  watts).  Similar  results  were  obtained  with  thicknesses 
greater  than  2  mm.  In  Table  2,  we  show-  the  power  of  the  laser  beam  incident  on  the 
layer  and  the  amount  of  radiation  reaching  the  detector  in  mW.  The  beam  suppression 
ratio  (BSR)  is  redefined  now  as  the  ratio  betw-een  output  power  divided  by  the  input 
power  and  is  given  in  Table  2.  The  BSR  ratio  is  the  reverse  of  the  Beam  Transmission 
Ratio  (BTR)  used  before.  We  chose  to  do  that  here  by  recognizing  it  is  easier  to  see 
how  much  energy  is  transmitted  through  the  layer  and  reaching  the  target.  It  is  clearly 
obvious  that  the  amount  of  energy  coming  out  of  the  layer  is  negligible  in  terms  of 
doing  any  harm  to  the  target,  i.e..  I0'6  to  I0'5  mW  (Tables  2  and  3). 
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TABLE  2 


BTR  of  NaCl  <  250^  Particle  Size,  2  mm  layer 

Time  of  exposure  of  the  layer  to  the  laser  beam  *  60  s 


Serial  # 

Power  Input 

Power  Output 

Beam  Transmission  Ratio 

mW 

mW 

BTR  ' 

1 

1000 

.0.008 

8.00E-06 

2 

1200 

0.008 

6.67E-06 

3 

1400 

0.009 

6.43E-06 

4 

1600 

0.010 

6.25E-06 

5 

0.C11 

6.11E-06 

6 

2000 

0.012 

6.00E-06 

7 

2400 

0.013 

5.42E-06 

8 

2800 

0.015 

5.36E-06 

9 

3200 

0.018 

5.63E-06 

10 

3600 

0.019 

5.28E-06 

1 1 

4000 

0.024 

6.00E-06 

12 

5600 

0.036 

6.43E-06 

1  3 

7600 

0.047 

6.18E-06 

14 

9600 

0.060 

6.25E-06 

1  5 

20000 

0.113 

5.65E-06 

I 
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TABLE  3 


BTR  of  NaCI  -  Particle  Size  =  25^111,  1  mm  Layer 
Time  of  exposure  of  the  layer  to  the  laser  beam  *  60  s 


Serial  # 

Power  Input 

Power  Output 

Beam  Transmission  Ratio 

mW 

mW 

BTR 

1 

1000 

0.027 

2.70E-05 

2 

1200 

0.034 

2.83E-05  1 

3 

2000 

0.056 

2.80E-05 

4 

4000  i 

0.114 

2.85E-05 

5 

7600  1 

0.221 

2.91  E-05 

6 

9600 

0.308 

3.21  E-05 

7 

20000 

0.836 

4.18E-05 

In  Table  4  we  show  the  results  for  layer  thickness  of  0.5  mm.  Again 
we  see  that  the  amount  of  energy  transmitted  through  the  layer  is  insignificant.  This 
fact  is  really  remarkable,  since  0.5  mm  is  equivalent  to  500Mm  thickness  of  a 
layer  of  approximately  25fim  particles,  which  in  turn  means  few  particles 
thick.  We  are  not  able  at  this  time  to  determine  values  for  layer  thicknesses  <  0.5  mm, 
for  the  difficulty  in  making  them.  We  would  like  to  determine  these  values  in  a 
continuation  study  that  would  enable  us  to  devise  methods  of  making  thin  (<  0.5  mm)  layers. 

TABLE  4 

BTR  of  NaCI  -  25jira  Particle  Size,  0.5  mm  layer. 

Time  of  exposure  of  the  layer  to  the  laser  beam  ■  60  5 


Serial  p 

Power 

Input 

Power  Oulput 

Beam  Transmission  Ratio 

mW 

mW 

BTR 

1 

1000 

0.130 

1.30E-04 

2 

1200 

0.164 

1.37E-04 

~3| 

1400 

0.185 

1.32E-04 

4  2000 

0.230 

1  15E-04 

5 

4000 

0.28< 

7.1 0E-05 

6 

5600 

0.280 

5.00E-05 

7 

7600 

0.311 

4.09E-05 

8 

9600 

0.361 

3.76E-05 

9 

15600 

0.996 

6.38E-05 

10 

17600 

1.830 

1.04E-04 

1 1 

19600 

3.360 

1 . 7 1 E- 04 

_L£j 

20000 

3.100 

1.55E-04 
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We  include  here  measurements  made  using  layers  of  particles  with  sizes  higher 
than  25  nm  in  order  to  see  if  larger  sizes  of  the  particles  will  affect  the 
reflectivity  of  the  layer  and  if  so  by  how  much: 

TABLE  5 

BTR  of  NaCI  -  40>iin  Particle  Size,  0.5  mm  Layer 
Time  of  exposure  of  the  layer  to  the  laser  beam  -  60  s 


Serial  # 

Power  Input 

1  Q 

Beam  Transmission 
Ratio 

mW 

mW 

BTR 

1 

1000 

0.16 

1 .60E-04 

2 

1200 

0.18 

1.50E-04 

3 

2000 

0.24 

1.20E-04 

4 

4000 

0.44 

1.10E-04 

5 

5600 

0.6 

1.07E-04 

6 

7600 

0.8 

1.05E-04 

7 

9600 

0.98 

1.02E-04 

8 

15600 

1.5 

9.62E-05 

9 

1  7600 

1.72 

9.77E-05 

10 

19600 

1.73 

8.83E-05 

1 1 

20000 

1.78 

8.90E-05 

TABLE  6 

BTR  of  NaCI  -  40jm  Particle  Sire,  0.5  mm  La>er 


Time  of  exposure  of  the  laser  to  the  laser  beam  -  60  s 


Serial  # 

Power 

Input 

Power  Output 

Beam  Transmission  Ratio 

mW 

mW 

BTR 

1 

1000 

0.32 

3  20E-04 

2 

1200 

3.00E-04 

3 

1400 

2.86E-04 

4 

2000 

0.49 

2.45E-04 

5 

4000 

0.95 

2.38E-04 

6 

5600 

1.32 

2.36E-04 

7 

7600 

1.74 

2.29E-04 

8 

9600 

1.93 

2.01E-04 

9 

2.91 

1.87E-04 

10 

mEBMi 

3  31 

1.68E-04 

1 1 

19600 

3  37 

1.72E-04 

1  2 

20000 

|  3  46 

1.73E-04 
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TABLE  7 


BTR  of  NaCJ  -  53/xm  Particle  Size,  0.5  mm  Layer 
Time  of  exposure  of  the  layer  to  the  laser  beam  *  60  s 


Serial  # 


Power  Input 


mW 


1000 


1200 


2000 


4000 


Power 

Output 


mW 


0.34 


0.42 


0.56 


1.07 


5600 


IE 


7600 


7  i  5600 


1  56C0 


17600 


10, 


19600 


111  20000 


1.33 


1.84 


2.18 


3.29 


3.64 


3.65 


3.68 


Beam  Transmission 
Ratio 


BTR 


3.40E-04 


3.50E-04 


2.60E-04 


2.68E-04 


2.38E-04 


2.42E-04 


2.27E-04 


2.11E-04 


2.0/E-04 


1.86E-04 


1.84E-04 


TABLE  8 

BTR  of  NaCI  -  90(xm  Particle  Size,  0.5  mm  Layer 
Time  cf  exposure  of  the  layer  to  the  laser  beam  -  60  s 


Se'.a  p 

Pc^er  irpu! 

Power  Output 

Beam  Transmission  Ratio 

mw 

mw 

BTR 

1 

1000 

0.32 

3.20E-04 

2 

1200 

0.36 

3.00E-04 

3 

2000 

0  62 

3.10E-04 

o 

o 

o 

'I 

*1 

1.15 

2.88E-04 

5j  5600 

1.62 

2.89E-04 

6 '  9600 

!  2  49 

2.59E-04 

7 |  15600 

3  64 

2.33E-04 

_ 1 

17600 

4  28 

9|  196CC 

4.34 

1 0f  20000  4  48 

2.24E-04 
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TABLE  9 


BTR  of  NaCl  -  125fiIH  Particle  Size,  0.5  mm  Layer 
Time  of  exposure  of  the  layer  to  the  laser  beam  ■  60  s 


Sena1  # 

Power  Input 

Power  Output 

Beam  Transmission  Ratiiy 

m  w 

m  w 

BTR 

1 

1000 

0.47 

4.70E-04 

2 

1  2C0 

0.5 

4.17E-04 

3 

20C0 

0.72 

3.60E-04 

4 

4000 

1.26 

3.15E-04 

5 

56C0 

1.68 

3.00E-04 

6 

96C0 

2.85 

2.97E-04 

7 

1  5600 

4.45 

2.85E-04 

6 

1  7600 

4.73 

2.69E-04 

9 

196  00 

4.92 

2.51E-04 

1C 

20C0C 

5.13 

2.57E-04 

TABLE  10 


BTR  of  NaCl  -  250/iIit  Particle  Siie,  0.5  mm  Layer 
Time  of  exposure  of  the  layt  \  .he  laser  beam  -  60  s 


Senal  #  Power  Input 


mw 


1000 


1200 


2000 


4000 


5600 


9600 


15600 


81  17600 


9 i  19600 


2CCOO 


Power 

Output 

Beam  Transmission 
Ratio 

mw 

BTR 

0.23 

2.30E-04 

0  33 

2.75E-04 

C.45 

2.25E-04 

0  91 

2.28E-04 

1.22 

2.18E-04 

1.96 

2.04E-04 

3.06 

1.96E-04 

3.19 

1.81E-04 

3.46 

1.77E-04 

3.48 

1.74E-04 

Part  III.  LASER/NUCLEAR  HARDENING  SCATTER  EFFICIENCY  MODEL 


This  Pari  presents  a  mathematical  model  for  the  properties  of  a  silica-based 
particle  layer  used  as  a  laser  shield.  Heat  flow  and  beam  attenuation  models  are 
derived,  then  combined  in  a  computer  program.  The  computer  model  predicts  maximum 
survivable  power  density  levels  for  varying  materials  and  shield  configurations.  In 
addition,  the  heat  flow  model  provides  temperature  profiles  across  the  shield  layer 
resulting  from  laser  energy  absorption,  and  the  beam  attenuation  model  provides  beam 
absorption,  reflection  and  transmission  as  a  function  of  layer  material  and  thickness. 

These  models  are  intended  to  provide  an  initial  theoretical  framework  to 
predict/extrapolate  performance  beyond  the  current  experimental  boundaries. 
Specifically,  the  model  can  be  used  to  extrapolate  the  performance  of  similar 
materials/coatings  in  the  mid-to-far  infrared  wavelengths,  and  to  predict  broadband 
performance.  Significant  work  remains  to  experimentally  verify  and  expand  the  results 
obtained  to  date,  and  to  further  develop  the  model  framework  begun  here. 

The  authors  of  this  model  continue  to  stress  the  wide  variety  of  applications  of 
this  laser  hardening  to  satellites,  satellite  solar  panels,  missile  boosters  and  nose 
cones.  SDI  deployed  systems,  and  others. 


MODEL  INPUTS 

1.  particle  size  (s). 

2.  particle  shape:  highly  irregular. 

3.  particle  refractive  index  (n). 

4.  particle  absorptivity  as  a  function  of  wavelength  a(  1). 

5.  powder  packing  fraction  (PF  -  %  of  solid  material). 

6.  particle  material  thermal-physical  properties: 

density  (r),  melting  temperature  Tm  .  thermal  diffusivitv  (  a)  which 
is  equal  to  thermal  conductivity  (k)  divided  by  the  heat  capacity 
<c).  and  density  (r). 

7.  layer  thickness  < x). 

8.  type  of  laser  beam  (CW  vs.  pulsed):  assumed  CW  for  duration  of  time  on  target. 

9.  wavelength  (I). 


10.  laser  power  density  (P) 


MODEL  OUTPUTS 


1.  amount  of  forward  scattering  or  transmission  of  laser  light 
through  the  layer. 

2.  amount  of  retroscattering  or  reflection  of  laser  light  from  the 
layer. 

3.  temperature  distribution  across  the  layer. 

4.  melting/damage  threshold  of  the  laser  shield. 

5.  bulk  absorption  through  the  layer. 


TEMPERATURE  PROFILE  MODEL 

Here  we  examine  the  flow  of  heat  in  the  coating  and,  in  particular,  the  resulting 
temperature  profile.  For  our  purpose,  which  is  to  obtain  a  worst-case  prediction,  a 
relatively  simple  mathematical  model  is  sufficient  to  describe  the  basic  features  of  the 
interaction  of  a  laser  beam  with  the  coating.  To  begin  with,  we  assume  that  the  coating  is 
an  infinitesimally  thin  homogeneous  and  isotropic  medium  with  thermal  conductivity  k, 

density  p,  and  specific  heat  c.  As  the  coating  is  essentially  two-dimensional,  the 
temperature  distribution  throughout  the  coating  is  governed  by  the  two-dimensional  heat 
equation 


IjL 

r  dr 


1  fT 

r2  M2 


+  ^g(r,8,l) 
k 


1 _dT_ 
a  dt 


(1) 


In  the  above,  T(r,Q,t)  is  the  temperature  at  the  point  (r,0)  (in  polar  coordinates)  at  lire  time 

t,  The  parameter  a  is  given  by  a  =  k/pc  .  The  function  g(r,0,t)  is  the  rate  of  heat 
generation  per  unit  volume. 

The  heat  equation  is  a  parabolic  differential  equation  and  needs  to  be  supplemented 
with  boundary  conditions  and  an  initial  condition,  Wc  assume  that  the  boundary  is  a  circle 
at  r=R  at  which  no  heat  can  enter  or  escape: 


dT_ 

dr 


r  =R 


=  0  . 


(2) 


At  t=0,  we  assume  that  Lite  layer  is  at  a  constant  temperature  T0,  Without  loss  of  generality, 
we  can  define  T  to  be  measured  relative  to  this  initial  temperature  (this  is  permissible  since 
both  the  heat  equation  (1)  and  our  boundary  condition  (2)  are  linear),  then  the  initial 
condition  is 
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T(r,e.0)-0.  (3) 

If  we  choose  our  coordinates  such  that  the  beam  is  centered  at  the  origin,  then 
azimuthal  symmetry  implies  that  T  is  independent  of  0  and  equation  (1)  reduces  to 


This  equation  is  solved  by  expanding  the  radial  dependence  of  T  (r.t)  in  terms  of  Bessel 
functions,  the  precise  form  of  the  expansion  is  controlled  by  the  boundary  condition  (2). 
Equation  (4)  reduces  to  a  linear,  first-order  ordinary  differential  equation  in  the  time  which 
is  easily  solved  in  conjunction  with  the  initial  condition  (3).  The  solution  takes  the  form 


T(r,t)  = 


m  *0  0 


where 


V2  •'„(&/) 

H  J0( A.A)  1 


and 


(5) 


(6) 


=  T"J  r)9(r-  n  • 

K  0  (7) 

The  functions  J0  are  zeroth-ordcr  Bessel  functions.  The  constants  are  the  roots  of  the 
equation 

dJAPr)  I 

*  _  n 


To  proceed  further,  we  need  to  choose  a  form  for  gfr.t).  Again,  for  our  purposes, 
the  simplest  choice  suffices.  We  model  the  laser  beam  as  a  point  source  of  constant 

integrated  power  y  located  at  r=0.  Mathematically,  this  is  accomplished  by  choosing  gfr,t) 
as 


(8) 


w  here  5u)  is  Dirac's  delta  distribution.  The  integral  over  r  in  (7)  can  then  be  performed: 


ay  1 
■s/ZnR  (An  R) 


(9) 
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This  then  allows  us  to  perform  the  integral  in  u  in  (5)  to  get 


T(r,t)  = 


-  e*-') 


(10) 


The  first  term  in  (10)  is  the  dominant  effect  of  the  boundary  and  represents  a  linear 
build  up  of  heat  in  lime.  The  remaining  infinite  series  is  the  exponential  approach  to  late- 
time  equilibrium  (obtained  by  setting  the  exponential  term  to  zero). 

We  evaluate  (10)  numerically  below,  but  a  useful  approximation  can  be  devised. 
Assume  the  medium  is  large  (PR»1),  and  that  we  observe  the  temperature  near  the  beam, 
i.e„  Pr«l,  then  using  the  asymptotic  forms  for  the  Bessel  functions  yields 


(ID 


The  cosine  term  in  ( 1 1 )  will  slowly  vary  between  ±  I .  The  sum  will  converge  due 

to  the  P term  which  can  be  approximated  by  (mit) w  ( In  (1 1)  we  renormalized  our 
distance  scale  such  that  R=l.)  For  small  r,  m  *R/4r,e.g.  For  R=1  cm  and  r=.01  cm, 

m(  =  25.  For  quartz,  a=.Q02  cm!/s,  so  for  t>  I  s,  pm^  >1000  will  give  an 
exponent  <  -2.  implying  that  the  exponential  term  can  be  neglected  for  m>10. 

For  higher  temperatures,  radiative  cooling  will  provide  a  small  but  significant 
contribution  to  the  thermal  distribution.  Inclusion  of  this  term  docs  not  readily  yield  an 
analytical  solution,  but  for  numerical  purposes  this  can  be  estimated  by  taking  the  time 
derivative  of  T(r,t),  and  subtracting  a  term  dependent  on  the  local  temperature.  In  general 

this  term  will  be  proportional  to  the  integral  of  the  emission  function  c(A)over  all 
wavelengths.  For  a  blackbody,  the  Planck  function,  and  the  integral  is: 


\b(X)dk  =  oT* , 

where  o  is  the  Stefan- Boltzmann  constant.  Thus  for  a  blackbody. 

dT  _  fry  y  '  a 
dt  ~  nk  +  2k  ^ e  ~  pc  1 

m  at 


Although  transparent  materials  are  poor  approximations  to  a  blackbody,  the  radiative  term 
is  small.  The  computer  model,  discussed  below,  ignores  the  radiative  cooling  term  in  its 
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present  implementation,  which  yields  a  more  conservative  estimate  of  shield  resistance  to 
melting. 

BEAM  ATTENUATION  MODEL 

This  section  deals  with  the  theory  of  the  relation  between  reflection,  transmission, 
and  absorption  to  layer  thickness.  We  will  define  the  parameters  upon  which  the  theory 
depends,  and  use  these  to  determine  a  thickness  which  will  give  a  reasonable  level  of 
protection  to  the  substrate. 

In  general,  w'e  wish  to  minimize  the  thickness  (x)  (and  thus  the  weight)  of  the 
shield,  while  maximizing  the  fraction  of  light  reflected  (R)  and  minimizing  the  fractions  of 
both  die  absorbed  (A)  and  transmitted  (T)  light.  Using  basic  principles,  we  can  derive  the 
equations  relating  these  quantities.  We  assume  that  for  some  thickness  x,  the  values  are 

known.  Then,  by  increasing  the  thickness  of  the  layer  by  an  additional  5x,  we  solve  for 

the  changes  in  these  quantities.  By  making  5x  small  enough,  w'e  can  assume  linearity  in 
reflection  and  absorption  so  that 

T  ( 5x )  =  1  -  cr<5x 
R{5x)  =  c^x 
A{5x)  -  Ca5x 

A  +  R  +  T  =  1  ,  nil 


where  ,  cT,  and  cA  arc  coefficients  of  reflection,  transmission,  and  absorption, 
respectively,  and  have  units  of  inverse  length.  These  are  related  such  that  cT  =ca  +ca  . 
Then,  using  multiple  scattering,  we  get 


T {x  +  Sx)  =  T(x)( 1-  cr6x)  +  T(x)R(x)cr6x( 1-  ct8x)  +  •  ■  ■ 
R{x+  5x)=  R(x)+  t\x)cr5x  +  t\x)  R{x)cr{5x)2  +  ■  •  • 


dT_ 

dx 


c,RT  -  cT  . 


(14) 


dR 

dx 


05) 


Equations  (13)  05)  can  be  solved  to  give: 


q(e2cx  - 1) 


(Q-1)(ecx-1) 
qecx  + 1 


where  c  =  (cTJ  -c%1  )'n  ,  and  q  =  (cA  +  c„  +  c)  /  ct  .  Define  x#  as  the  thickness  where  half 
the  light  is  transmitted^  then  x0  =  1  when  c,  =  1.  If  there  were  no  absorption,  then 


In  the  limit  of  an  infinitely  thick  layer, 

A  =  ( 2A ^  (20) 

where  the  subscript  "o"  denotes  the  absorption  of  a  layer  of  thickness  x0 .  A0  is 
approximately  equal  to  ca  xowhcn  ca  is  small. 

Because  absorption  increases  with  increasing  thickness  (to  a  limit)  and  transmission 
decreases,  some  reasonable  criterion  must  be  set.  To  minimize  the  sum  of  absorption  plus 
transmission  seems  reasonable  at  first,  but  this  sum  will  always  decrease,  albeit 
insignificantly,  as  the  thickness  increases.  Let  us  define  the  "ideal"  thickness  (x.),  as  that 
where  A  =  T.  Beyond  this  thickness,  the  absorbed  energy  dominates  the  damage  potential. 
This  is  not  to  say  that  increasing  the  thickness  will  increase  the  damage  potential,  since  the 
absorbed  energy  per  unit  mass  will  decrease.  Nonetheless,  setting  T  =  A  will  give  us  a 
suitable  criterion  to  work  with.  Substituting  the  expressions  in  equations  (4)  and  (6)  and 
solving  for  x  gives: 


=  7-Wta  +  i)  +  V(q  +  I)2  -  1J-  I  n  g  } 


Given  ca  /cT ,  this  can  be  solved  for  x.  A  simple  approximation  is  also  available. 
Absorption  is  ncarl)  linear  with  thickness  for  x  <  xo  and  transmission  approximately 
inversely  proportional  to  thickness.  So  we  can  get  a  quick  estimate  of  the  thickness  for 
maximum  beam  suppression: 


*  =  *cAi 


For  example,  if  Ao  =  10"  .  then  the  estimated  ideal  thickness  is  100  x0 ,  at  which 
point  both  absorption  and  transmission  are  1%  of  the  total  radiation.  The  numbers  for  the 
exact  solution  come  out  slightly  better,  absorption  =  transmission  =0.81%  for  a  thickness 
of  92.6  x#. 

We  can  minimize  the  thickness  of  the  layer  by  minimizing  xfl.  This  can  be  done  to  a 
limited  extent  by  scaling  down  the  sizes  of  the  panicles.  Eventually,  in  the  limit,  the 
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particles  become  smaller  than  the  wavelength  of  the  radiation,  and  their  scattering  properties 
change.  Still,  if  x0  is  10  panicle  thicknesses,  then  use  of  1  |im  particles  gives  a  total  layer 
thickness  of  less  titan  1  mm  in  the  above  example. 

There  is  a  theoretical  reason  to  believe  that  when  particle  size  is  decreased,  the 
corresponding  decrease  in  x0  is  better  titan  linear.  This  is  due  to  the  fact  that  the  smaller  the 
particle,  llie  less  peaked  are  the  scattering  curves  in  the  forward  direction.  It  would  take  a 
several-uarticleS'thick  layer  of  large  particles  to  generate  the  equivalent  backscattcr  of  a 
monolayer  of  small  particles. 

MODEL  IMPLEMENTATION 

The  model  equations  derived  above  have  been  incorporated  into  a  computer 
program.  Fortran  source  code  for  this  program  is  presented  in  Appendix  C.  Most  input 
parameters  arc  requested  at  program  initiation.  However,  the  thermal-physical  constants  of 
the  material  must  be  changed  in  the  code  itself,  if  a  material  other  than  silica  is  used.  See 
A  ppendix  A  for  the  parameter  input  sequence  and  default  values,  The  program  iterates  to 
compute  the  power  density  required  to  bring  the  layer  to  99.5%  of  the  specified  melting 
temperature,  w  ithin  the  specified  time  the  beam  is  on  target.  This  power  density  is  then 
presented  as  die  resulting  "maximum  power  density"  die  shield  can  withstand  given  die 
specified  input  conditions.  See  Appendix  B  for  a  sample  output  of  the  model. 

The  effect  of  certain  properties  of  the  material  on  the  transmission  and  heating  of  the 
layer  are  presently  unknown.  Many  of  these  relationships  arc  planned  subjects  of  our 
proposed  experimental  study.  To  facilitate  the  creation  of  a  model,  we  have  estimated  the 
functional  dependences  described  below. 

The  effect  of  panicle  size  was  decoupled  into  two  separate  forms.  The  geometric 
assumption  was  that  an  n-panicle  thick  layer  scatters  with  equal  efficiency,  independendy 
of  the  size  of  the  panicle.  Superimposed  on  this  was  the  assumption  that  the  closer  the  size 
of  die  panicle  is  to  die  wavelength  of  the  incident  beam,  the  more  efficiendy  tt  scatters  away 
from  tiie  forward  direction  cf.  van  de  Hulst,  1957^.  The  difference  includes  both  an 
increase  in  die  ratio  of  true  scattering  cross  section  to  geometnc  scattering  cross  section  (a 

factor  ol  two  at  s=k;Spitzer,  I968',)tmd  an  effect  due  to  the  increase  in  isotropy  of  the 
scattering  curve  of  small  particlcsover  large  particles.  This  was  approximated  by 


cr  -  j-(y  +logy)’1 

(24) 

where  Ct  ;s  inversely  proportional  to  the  thickness  needed  for  5U%  reflection,  so  that  a 
larger  value  of  Cj  means  a  more  efficient  scattcrer. 

The  index  of  refraction  was  also  expected  to  play  a  role  in  die  scattering  efficiency, 
based  on  Mic  scattering  data.  This  was  modeled  by 


n  -  1 
n 


to  give  a  number  varying  between  zero  (n=l)  and  one  (n  infinite). 


(25) 
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Finally,  n  was  assumed  that  the  maximum  density  of  interfaces  is  provided  by  a 
layer  of  75%  solid  material,  with  a  higher  density  allowing  faces  to  overlap  and  a  lower 
density  allowing  excess  air  (or  vacuum)  space.  c-y  must  go  to  zero  when  the  layer  is  all  air 
(PF=0)  or  solid  material  (?F=1),  Using  a  cubic  to  model  this  property  and  combining 
these  expressions  gives: 


c  _  9P£_{n  -  1)(  -  8PF2  -t-  1 1PF  -  3)  (25) 

r  n  (A  +  s  log-j) 

where  the  9  is  a  normalizing  term,  based  on  a  conservative  fit  to  our  previous  experimental 
results. 


The  program  will  allow  either  the  thickness  or  the  maximum  percentage 
transmission  to  be  input  If  both  arc  set  to  zero,  the  code  will  calculate  the  default  thickness 
where  Transmission  =  Absorption. 

Radiative  cooling  has  been  eliminated  in  the  cument  implementation  by  setting  the 
Stefan-Boltzmann  constant  to  zero.  This  gives  a  lower  limit  on  the  resulting  maximum 
survivable  power  density. 

SAMPLE  PROGRAM  RESULTS 

We  like  to  acknowledge  here  the  efforts  of  Dr.  Edwin  T.  Rusk  and  Dr.  Charles  G. 
Tome  at  FIT  in  developing  the  mathematical  treatment  and  the  computer  coding  of  this 
model.  Appendix  B  provides  a  sample  output  of  computer  model  results.  The  example 
temperature  oistriDuuon  is  foi  2  seconds  in  steps  of  .01  s  It  was  generated  for  a  beam  of 
.01  cm  diameter  centered  on  a  layer  1  cm  in  radius  and  set  the  maximum  temperature  at  the 
melung  point  of  quartz.  Otherwise,  the  default  parameters  of  the  model  were  used.  These 

include  an  absorption  coefficient  of  SxlO'Vcm,  panicles  1pm  in  djameter,  a  wavelength  of 
488  n:n,  an  index  of  refraction  of  1 .46,  and  a  packing  fraction  of  75%.  See  Figure  9  for  a 
graphic  presentation  of  the  results. 

The  program  iterates  on  beam  power  density'  to  produce  a  temperature  of  1873.6K 
(99.5%  of  1 883K,  the  melting  point)  after  2  seconds.  Thus,  the  computed  maximum 
power  density  of  34.6  MW/cm-  is  the  power  density  the  shield  can  withstand  for  two 
seconds  without  melting.  Note  that  the  ideal  layer  thickness  computation  to 
provide  Transmission  =  Absorption  resulted  in  a  layer  of  0.932mm  and  a 
transmission/absorption  of  0.04%. 
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Figure  9;  plot  of  temperature  (Kelvin)  at  beam  edge  (.005  cm 
radius)  over  the  2  second  beam  exposure,  as  presented  in  Attachment  B. 
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CONCLUSION 


These  experiments  proved  that  it  is  possible  to  reject  more  than  99.9%  of  the 
argon  laser  beam  using  highly  pure  natural  occurring  beach  sand  in  Florida.  This 
rejection  was  accomplished  by  packing  highly  irregular  particles  in  thin  layers  (0.5, 

1.0,  1.5,  and  2.0  mm),  which  then  creates  a  multiple  scattering  medium  for  the  laser 
beam.  The  absorption  coefficient  of  these  highly  pure  silica  particles  is  of  the  order 
of  4  x  10”^  cm  which  makes  them  resistant  to  melting  by  the  laser  energy.  This 
rejection  of  the  radiation  away  from  the  target  is  feasible  between  0.22  to 
2.4  /im  wavelength.  In  this  research  task,  we  embarked  on  extending  this 
method  to  reject  more  than  99.99%  of  the  radiation  from  the  beam  of  a  20  Watt  CO 2  CW 
laser,  using  NaCl  (natural  occurring  salt).  This  extends  our  methodology  to  reject 
radiation  away  from  the  target  from  0.22  to  15/im  region  of  the  spectrum. 

The  CW,  CO2  laser’s  beam  transmission  ratio,  BTR  (the  ratio  of  the  energy 
transmitted  to  the  target  divided  by  that  reflected  away  from  the  target,  see  Figure  1), 
is  the  reverse  of  BSR  that  we  used  in  the  initial  experiments  with  silica  particles.  We 
preferred  to  use  BTR  since  it  shows  readily  how  small  is  the  energy  reaching  the 
target.  Tables  2  to  10  show  clearly  the  insignificant  amount  of  energy  reaching  the 
target.  We  were  pleasantly  surprised  as  to  how  well  NaCl  particles  reflected  the  energy 
from  the  CO2  laser  using  very  thin  layers  <  0.5  mm,  compared  to  the  silica  particles 
layers.  We  suspect  that  the  reason  for  that  is  the  long  wavelength  (10.6  fim  ) 
of  the  photons  coming  out  of  the  C02  laser,  and  their  increased  inability  to  emerge 
out  of  the  layer  compared  to  the  photons  in  the  visible. 

Natural  NaCl  has  an  absorption  coefficient  of  approximately  10  cm  '.  This 
very  low  absorption,  from  2  to  15)jro  wavelength,  ensures  that  this  laser 
shield  will  not  melt  or  sublimate  easily.  We  have  subjected  layers  of  thicknesses  0.5, 

1.0,  1.5,  2.0,  and  4.0  mm  to  energy  density  of  approximately  1.5  KW/cm2  (the  limit 
of  our  focused  CW  CO2  laser  beam)  and  did  not  notice  any  melting.  FIT  likes  to  test 
for  melting  at  higher  energy  densities  and  will  be  awaiting  the  use  of  the  2.5  KW,  C02 
laser  at  Natick  RD  &  E  Center  for  such  tests. 

We  also  experimented  with  adhering  the  highly  irregular  NaCl  particles  in  the 
layer,  using  commercial  adhesives.  The  best  adhesive  we  have  used  so  far  is  the 
commercially  available  polyurethane  spray  adhesive.  The  application  of  our 
laser/nuclear  shielding  is  fairly  easy  for  rigid  wall  hardening.  The  application  of  our 
shield  on  body  armor  is  more  complicated  due  to  the  fact  that  it  must  be  applied  in  a 
way  that  keeps  the  body  armor  flexible  and  comfortable.  Using  this  coating  material  on 
body  armor  is  subject  to  more  research  if  the  needs  of  the  Army  demand  hardening 
against  CO->  laser  threats  and  nuclear  threats  as  well. 

Currently,  we  are  exploring  the  use  of  our  coaling  to  reduce  the  solar  load  off 
Army  containers  and  canisters  (rigid  walls).  The  same  method  with  which  we  reject  the 
radiation  from  laser  beams  can  be  applied  to  reflecting  most  of  the  solar  spectrum  thus 
making  the  environment  cooler  in  hot  battlefields  such  as  those  encountered  in  the 
Desert  Shield  and  Desert  Storm  operations  in  Saudi  Arabia  and  southern  Iraq.  NaCl  will 
reflect  the  infrared  part  of  the  solar  spectrum  and  natural  highly  pure  silica  (Florida 
beach  sand)  rejects  the  solar  spectrum  from  0.22  to  2.4fim.  We  have  not  tested 
this  coating  against  the  solar  spectrum  below  0.22  fxm;  however,  application  of 
this  coating  on  rigid  walls  and  exposing  it  to  sunlight  should  demonstrate  its 
effectiveness  in  this  region  of  the  solar  spectrum. 
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Polyurethane  as  adhesive  proved  to  be  flexible  on  metal  surfaces.  It  resists 
cracking  due  to  twisting^ or  bending  of  the  metal.  The  coating  is  light  weight, 
approximately  1.2  kg/m2  when  0.5  mm  in  thickness. 

We  hope  that  our  findings  from  this  research  task  will  prove  to  be  very  useful 
for  the  U.S.  Army's  defensive  applications  in  the  battlefield  theater  against  laser  and 
nuclear  threats. 
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Appendix  A:  Parameter  Input  Sequence 


ENTER  THE  MAXIMUM  ALLOWABLE  TEMPERATURE ,  ELSE  2200. 00K 
1883  . 

THE  DIAMETER  OF  THE  BEAM ,  ELSE  0.100  cm 
.01 

THE  DIAMETER  OF  THE  LAYER,  ELSE  1.00  cm 

2. 

THE  TIME  THE  BEAM  IS  ON  TARGET,  ELSE  10.00  sec 

2. 

THE  TIME  INCREMENT,  ELSE  0.1000  sec 

.01 

THE  ABSORPTION  COEFFICIENT,  (ELSE  0.5E-03  per  cm) 

THE  LAYER  THICKNESS,  (ELSE  " IDEAL” ) 

THE  PARTICLE  SIZE,  (ELSE  1.00  um) 

THE  WAVELENGTH,  (ELSE  0.4880  um) 

THE  l -AGE  OF  SOLID  MATERIAL,  (ELSE  75.  %) 

THE  INDEX  OF  REFRACTION,  (ELSE  1.460) 

THE  MAXIMUM  ALLOWABLE  TRANSMISSION  %age,  IF  ANY. 


3i 


Appendix  B:  Sample  Model  Output 


?  OR  Xi=  0 . 9  3  2  nun ,  0.04006  t  IS  TRANSMITTED,  AMD  0.04006  %  IS  ABSORBED 


TEMPERATURE  PROFILE 


RADI  I  ( cm .  ) 

l  L  '  1 


TIME  (s. 

)  0.00 

0 . 14 

0 . 17 

0.20 

0.25 

0.33 

0.50 

1.00 

0.010 

483  .  1 

300.1 

300.  1 

300.  1 

300.1 

299  .  9 

300  .  1 

300.0 

0 . 020 

601 . 1 

300.1 

300.  1 

300 . 1 

300  .  1 

299.9 

300.1 

300.0 

.  0.030 

686.6 

300.  1 

300 . 2 

300 . 2 

300 . 1 

299.8 

300.2 

300.0 

0.040 

753  .  5 

300 . 1 

300.2 

300.2 

'300.1 

299.8 

300.3 

300.0 

0.050 

808 . 4 

300  .  1 

300.2 

300.2 

300 . 1 

299  .  8 

300.3 

300.0 

0.060 

854.9 

300  .  1 

300.2 

300 . 2 

300.  1 

299.8 

300.3 

300.0 

0.070 

895 . 3 

300.  1 

300.2 

300.2 

300  .  1 

299.8 

300.4 

300.0 

0.080 

930.9 

300.1 

300.2 

300 . 2 

300  .  1 

299.8 

300.4 

300.0 

0.090 

962.8 

300.  1 

300.2 

300.2 

300.1 

299  .  8 

300.4 

300.0 

0.100 

991.7 

300 . 1 

300.2 

300.2 

300 . 1 

299.8 

300.4 

300.0 

0 .120 

1042.3 

300.1 

300.2 

300.2 

300  .  1 

299.8 

300.4 

300.0 

0.160 

1123  .  7 

300.1 

300.2 

300.2 

300.1 

299.8 

300.5 

300.0 

0.200 

1187.8 

300.1 

300.2 

300.2 

300 . 1 

299.8 

300 . 5 

300.0 

0.240 

1240.7 

300.  1 

300 . 2 

300 . 2 

300  .  1 

299.8 

300 . 5 

300.0 

0 .280 

1285.7 

300.1 

300.2 

300 . 2 

300.1 

299.8 

300.5 

300.0 

0.320 

1325.0 

300.1 

300.2 

300.2 

300.1 

299.8 

300.5 

300.0 

0 . 360 

1359.7 

300.1 

300 . 2 

300 . 2 

300.  1 

299.8 

300.5 

300.0 

0.400 

1390 . 9 

300.  1 

300.2 

300.2 

300  .  1 

299  .  8 

300.5 

300.0 

0 .440 

1419.1 

300.1 

300.2 

300.2 

300.  1 

299.8 

300.5 

300.0 

0.480 

1445.0 

300.1 

300.2 

300.2 

300.1 

299  .  8 

300.5 

300.0 

0.520 

1468.8 

300  .  1 

300.2 

300.2 

300. 1 

299  .  8 

300 . 5 

300.0 

0.  560 

1490.9 

300.1 

300.2 

300 . 2 

300 . 1 

299  .  8 

300.5 

300.0 

0.600 

1511.5 

300 . 1 

300 . 2 

300 . 2 

300. 1 

299.8 

300.5 

300.0 

0.640 

1530.8 

300.2 

300.2 

300.2 

300.1 

299  .  8 

300.5 

300.0 

0.680 

1548.9 

300.2 

300.2 

300.2 

300.  1 

299 . 8 

300.5 

300.0 

0.720 

1566.0 

300.2 

300.2 

300 . 2 

300.  1 

299 . 8 

300.5 

300.0 

0.760 

1582 . 2 

300.3 

300.2 

300.2 

300. 1 

299  .  8 

300.5 

300.0 

0.800 

1597.6 

300.4 

300.2 

3  00.2 

300  .  1 

299 . 8 

300.5 

300.0 

.0.840 

1612.2 

300.5 

300.2 

300.2 

300 . 1 

299.8 

300.5 

300 . 0 

0 .880 

1626.2 

300 . 6 

300 . 2 

300.2 

300.1 

299.8 

300.5 

300.0 

0.920 

1639.5 

300.7 

300.3 

300.2 

300.1 

299.8 

300.5 

300.0 

0.960 

1652.3 

300.9 

300.3 

300.2 

300.1 

299 . 8 

300.5 

300.0 

1.000 

1664  .  6 

301.1 

300.3 

300.2 

300. 1 

299  .  8 

300.5 

300.0 

1.040 

1676.4 

301.3 

300 . 4 

300.2 

300.  1 

299.8 

300.5 

300.0 

1.080 

1687.7 

301.5 

300 . 5 

300 . 2 

300.  1 

299 . 8 

300.5 

300.0 

1.120 

1698.7 

301.8 

300.5 

300.2 

300.1 

299 . 8 

300.5 

300.0 

1 .160 

1709.3 

302  .  1 

300.6 

■  300.2 

300 . 1 

299.8 

300.5 

300.0 

1.200 

1719.5 

302.4 

300 . 7 

300 . 2 

300 . 1 

299 . 8 

300 . 5 

300.0 

1.240 

1729 . 3 

302.7 

300.8 

300.2 

300.  1 

299  .  8 

300.5 

300.0 

1.280 

1738  .  9 

303 . 1 

300.9 

300.2 

300 . 1 

299  .  8 

300.5 

300.0 

1.320 

1748.2 

303.5 

301.0 

300.3 

300  .  1 

299 . 8 

300.5 

300.0 

1.360 

1757.2 

303.9 

301.2 

300.3 

300.1 

299  .  8 

300.5 

•  300.0 
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1.400 

1765.9 

304  .  4 

301.3 

300 . 3 

300.  1 

299  .  8 

300.5 

300 . 0 

1 .440 

1774  .  4 

304 . 8 

301.5 

300  .  3 

300  .  1 

299  .  8 

300.5 

300.0 

1.480 

1782.7 

305 . 3 

301.7 

300.4 

300.  1 

299  .  8 

300.5 

300.0 

1 .520 

1790.7 

305.9 

301.9 

300.4 

3  00.1 

299.8 

300 . 5 

300.0 

1.560 

1798.6 

306.4 

302  .  1 

300.4 

300  .  1 

299  .  0 

300.5 

300  .  0 

1  .  600 

1806  .  2 

307  .  0 

302  .  3 

300 . 5 

300  .  1 

299.8 

300 . 5 

300.0 

1  .  640 

1813.7 

307.6 

302  .  5 

300 . 5 

300  .  1 

2  9  9  .  8 

300  .  5 

300.0 

1  .  680  ' 

1820 . 9 

308 . 2 

302  .  8 

300  .  6 

300  .  1 

2  99.8 

300  .  5 

300.0 

1.720 

1828.0 

308  .  8 

303  .  0 

300  ..7 

300  .  1 

2  9  9  .  8 

300  .  5 

300.0 

1  .760 

1  8  35.0 

309 . 5 

303.3 

300.7 

300 . 2 

299.8 

300  .  5 

300.0 

1  .  800 

1841.8 

310.1 

303  .6 

300.8 

300  .  2 

299.0 

300 . 5 

300. 0 

1.840 

184  8.4 

310.8 

303 . 9 

300  .  9 

300  .  2 

299  .  8 

300 . 5 

300 . 0 

1.880 

1  854.9 

311.5 

304  .  2 

301.0 

300  .  2 

299  .  8 

300 . 5 

300.0 

1.920 

186  1.3 

312.3 

304  .  6 

30  1 .  1 

300 . 2 

299  .  0 

300  .  5 

300.0 

1 .960 

1067 . 5 

313.0 

304 . 9 

301.2 

300 . 2 

2  99.8 

300 . 5 

300.0 

2  .  000 

1873.6 

313.8 

305.3 

301.3 

300  .  2 

299  .  8 

300.5 

300.0 

THE  CENTRAL  TEMPERATURE  IS:  1873.59  K 

THE  MAXIMUM  POWER  DENSITY  IS  :  28.8222  MW/ cm2 

IT  TOOK  '  3  STEPS  TO  FIND  THE  POWER 


r' 
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Appendix  C:  FORTRAN  Source  Code 


THIS  PROGRAM  GENERATES  A  MODEL  Or  THE  SCATTERING  PROPERTIES  Or  A  LAYER 
PROGRAM  LAYER 

IMPLICIT  REAL* 8  ( A-H , J , L , N , 0-Z ) 

REAL*  8  XR( 10) , HR ( 10 ), H( 10000 ), K,T{ 10, 1000 ) ,TIME( 1000  ) 

LOGICAL  NOTMAX ,  NOX ,  HIGH,  POWER 
EXTERNAL  JO 
COMMON'/  BETA  /  B 

Z  RO=  0 . DO 
HAF=  . 5D0 
ONE-  1 . DO 
TWO=  2 . DO 

PI=  3. 14 1592653509793 DO 

THESE  PARAMETERS  CAN  BE  VARIED,  DEPENDING  ON  THE  MATERIAL 

OPEN ( 10, STATUS' 'NEW' , F I LE  = ' LAYER . DAT ' ) 

K=  . 26D0 
SIGMA®  5 . 6  7  D—  8 

SIGMA*  ZRO 

1134=  8  .  1D9 
1134=  ZRO- 
RHO=  2.65D3 
CTH=  794. 

ALPHA*  K  /  (RHO*CTH) 

113=  3.D2 

I  RM=  8 

IRM0=  IRM 

IR2=  NINT(  4.5  *  IRM  ) 

CALL  BET A0 

CA0=  5.D-4 

S0=  ONE  , 

L0=  . 4  8  8  DO 

PP0=  75. DO 
NO*  1.46D0 
TMELT0*  2200. 

R0=  ONE 
RB0  =  . 1D0 
TYM0*  1 . D 1 

DT0=  l.D-1 
X0=  ONE 
Pl=  1.D2 

WRITE (6,20)  TMELT0 
READ ( 5 , 120)  TMELT 

I F (  TMELT  .LE.  ZRO  )  TMELT*  TMELTO 
TMELT0*  TMELT 


WRITE (6, 32)  RB0 
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READ (5,120)  RB 
I F (  RB  .LE.  ZRO  )  RB  =  RBO 
RBO=  RB 

RB=  RB  *  . 5D-2 

C 

WRITE (6,21)  R0 
READ ( 5 , 120)  R 
I F (  R  .LE.  ZRO  )  R=  RO 
RO=  R 

R=  R  *  .  5D-2 

C 

WRITE (6,22)  TYMO 
READ (5,120)  TYM 
I F {  TYM  .LE.  ZRO  )  TYM=  TYMO 
TYMO=  TYM 

C 

WRITE (6,29)  DTO 
READ ( 5 , 120)  DT 
I F (  DT  .LE.  ZRO  )  DT=  DTO 
DT0=  DT 

C 

WRITE (6,51)  CAO 
READ ( 5 , 120)  CA 
I F {  CA  .LE.  0 . DO  )  CA=  CAO 
CAO  =  CA 


WRITE (6,92) 

READ (5,120)  X 
NOX=  X  .LE.  0 . DO 
X=  X  *  .  1 

C 

WRITE ( 6 , 9  3 )  SO 

READ (5, 120)  S 

I F (  S  .LE.  0 . DO  )  S=  SO 

.  so-  s 

S=  S  *  l.D-4 

C 

WRITE (6,94)  LO 

READ ( 5 , 120)  L 

IF (  L  .LE.  0 . DO  )  L=  LO 

L0=  L 

L=  L  *  l.D-4 

C 

WRITE ( 6 , 9 5  )  PPO 
READ (5,120)  PP 
I F {  PP  .LE.  0 . DO  )  PP=  PPO 
PPO=  PP 

PP=  PP  *  l.D-2 

C 

WRITE(6,96)  NO 
READ (5, 120)  N 
I F (  N  .LE.  0 . DO  )  N=  NO 
NO-  N  ' 


C 
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WRITE (6,97) 

READ( 5,120)  TMAX 
NOTMAX=  TMAX  . LE .  0.D0 
TMAX=  TMAX  /  1.D2 

Y=  9. DO  *  PP  *  (  1 1 . DO  *  PP  -  8 . DO  *  PP*  PP  -  3 . DO ) 

Y=  (N-ONE)  *  Y  /  (L/S  +  DLOG10(S/L)) 

CT=  Y  /  (  N  *  S  ) 

CA=  1.D1  *  CA 
CR=  CT  -  CA 

C=  DSQRT(  CA  *  (CT+CR)  ) 

Q«  (CA  +  CR  +  C)  /  CR 
QPl=  Q  +  1 
QM1-  Q  -  1 
Q2M=  QP1  *  QM1 

XI=  QP1  +  DSQRT(  QP 1 *QP 1  -  1  ) 

XI=  DLOG (  XI  )  -  DLOG(  Q  ) 

X I =  XI  /  c 

EX=  DEXP(  C  *  XI  ) 

TI=  1.D2  *  Q2M  *  EX  /  (  Q  *  Q  *  EX  *  EX  -  1  ) 

AI=  1.D2  *  QM1  *  (  EX  -  1  )  /  (  Q  *  EX  +  1  ) 

XI=  1.D1  *  XI 

X0=  XI 

ABS=  AI  *  l.D-2 
WRITE( 6,96)  XI,  TI,  AI 
WRITE ( 10,98)  XI,  TI,  AI 

I F (  NOX  )  GO  TO  5 
EX=  DEXP(  C  *  X  ) 

AX*  1.D2  *  QM1  *  (  EX  -  1  )  /  (  Q  *  EX  ♦  1  ) 

TX*  1.D2  *  Q2M  *  EX  /  (  Q  *  Q  *  EX  *  EX  -  1  ) 

X=  X  *  1.D1 
X0=  X 

ABS=  AX  *  l.D-2 
WRITE( 6,99)  X,  TX,  AX 
WRITE( 10,99)  X,  TX,  AX 

5  IF (  NOTMAX  )  GO  TO  10 

XT=  Q2M*Q2M  +  4 *Q*Q*TMAX*TMAX 
XT*  Q2M  +  DSQRT(  XT  ) 

XT=  XT  /  (  2  *  Q  *  Q  *  TMAX  ) 

XT=  DLOG (  XT  )  /  C 
EX=  DEXP (  C  *  XT  ) 

AT=  1.D2  *  QM1  *(EX-1)/(Q*EX+1) 

TT2  =  1.D2  *  Q2M  *  EX  /  (  Q  *  Q  *  EX  *  EX  -  1  ) 
XT=  XT  *  i.Dl 
X0-  XT 

ABS*  AT  *  l.D-2 

WRITE (6,99)  XT,  TT2 ,  AT 

WRITE (10,99)  XT,  TT2 ,  AT 

3t 


;  THIS  SECTION  CALCULATES  THE  MAXIMUM  BEAM  POWER  BEFORE  MELTING 
C 

10  X-  X0  *  1 .D-3 
RHO=  RHO  *  PP 
BF=  RB  /  R 

P2=  PI  /  (PI  *  RB  *  RB) 

C 

DO  45  1=2, I RM 

XR ( I ) =  ONE  /  ( ONE - I RM- I ) 

45  RR( I ) «  XR( I )  *  R  *  1 . D2 

XR (  1 )  =  XR(IRM)  *  BF 

RR ( 1 ) =  XR( 1)  *  R  •  1.D2 

C 

ITM=  NINT(  TYM  /  DT  ) 

RAD1  =  SIGMA  *  DT  /  (  RHO  *  CTH  *  X  ) 

EXP 1 =  ALPHA  *  DT  /  (  R  *  R  ) 

TEM0=  ALPHA  *  DT  *  BF  *  BF  *  ABS  /  (K  *  X) 

C 

C  LOOP  OVER  POWER,  TO  FIT  FINAL  TEMPERATURE 
C 

Tl=  ZRO 
TG=  1.D7 
PG=  1.D13 
POWER=  .FALSE. 

HIGH=  .FALSE. 

IRM=  1 
C 

DO  350  I P=  1 , 9 9 9 
TEM1 =  TEMO  *  P2 

C 

DO  50  IT= 1 , ITM 
DO  50  IR= 1 , I RM 
50  T ( IT, I R ) =  H3 
C 

C  LOOP  OVER  RADIUS 
C 

IF (  POWER  )  I RM  =  IRM0 
DO  300  I R= 1 , IRM 
C 

C  LOOP  OVER  TIME 
C 

TOLD=  H3 

DO  300  IT= 1 , ITM 
TIME (  IT )  =  IT  *  DT 
BEXP  =  EXP  1  *  IT 

RAD=  ZRO 
C 

C  THIS  SECTION  CALCULATES  THE  SUMMATION  TERM 
C 

100  SUM=  ZRO 

DO  200  1=1,9999 

DSUM=  JO ( I , XR ( I R ) )  *  DEXP (  -BEXP * B ( I ) * B ( I )  ) 

SUM  =  SUM  +  DSUM 

IF (  DABS (DSUM)  . LT .  DABS ( 1 . D- 6  *  SUM )  )  GO  TO  2  50 
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n  n  n  n  n  n  n  n  n  non 


200  CONTINUE 

WRITE (6,28)  IP 

THIS  SECTION  CALCULATES  THE  RADIATIVE  LOSSES 

250  TEMP=  TEM1  *  (  ONE  +  SUM  ) 

TT=  TOLD  +  TEMP  *  HAF 

IF (  TT  .GT.  1.D8  )  WRITE(6,*)  IR,  IT,  TOLD,  TEMP 

T4=  TT**4 

RAD=  RADI  *  (T4  -  H34  ) 

I F (  RAD  .GT.  TT  )  RAD=  TT 
TOLD=  TEMP  -  RAD  +  TOLD 

T( IR, IT) =  TOLD 

300  CONTINUE 

P3=  P2  *  l.D-10 

WRITE ( 6 , * )  IP,  P3,  RAD 

WRITE ( 6 , * )  (T(IR,1),IR=1, IRM) 

WRITE ( 6 , * )  (T( IR, ITM) , IR= I , IRM) 

T2=  T ( 1 , ITM) 

IF (  POWER  )  GO  TO  375 
IF (  IP  .GT.  1  )  GO  TO  325 
TL=  T( 1 , ITM) 

PL-  P2 

P2=  5. DO  *  (TMELT  -  H3)  *  P2  /  (TL  -  H3) 

GO  TO  350 

TEST  FOR  TEMPERATURE  WITHIN  CRITICAL  VALUE 

325  POWER =  T2.LE. TMELT  .AND.  T2 .GT. TMELT* . 99 

IF (  POWER  )  GO  TO  350 
IF (  HIGH  .OR.  T2.GT. TMELT  )  GO  TO  330 
P2=  1.D1  *  P2 

GO  TO  350 

33C  HIGH=  .TRUE. 

IF (  T2  .GT.  TMELT  )  TG=  T2 

IF(  T2  -GT.  TMELT  )  PG=  P2 

IF(  T2  .LT.  TMELT  )  TL=  T2 

IF (  T2  .LT.  TMELT  )  PL=  P2 

DTEM=  TG  -  TL 

DTMELT=  TMELT*. 995  -  TL 
DP=  PG  -  PL 

P2*  PL  ♦  DP  *  DTMELT  /  DTEM 
350  CONTINUE 

WRITE  TEMPERATURE  DISTRIBUTION  AND  MAXIMUM  POWER 
375  CONTINUE 

WRITE ( 10,26)  ( RR ( I ) ,1=1, IRM) 

DO  400  IT=1,  10 

IF (  IT  .GT.  ITM  )  GO  TO  500 
400  WRITE( 10,121)  TIME(IT),  (T(IR,IT),  IR=1,IRM) 
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IDT=  MAX (  ITM/50  ,  1  ) 

DO  450  IT=11,ITM 

IF (  IT  .GT.  ITM  )  GO  TO  500 

IF (  MOD( IT, IDT)  . NE .  0  )  GO  TO  450 

WRITE( 10/ 121)  TIME ( IT) ,  (T(IR,IT),  IR=1,IRM) 

450  CONTINUE 

C 

500  WRITE( 6,25)  T2,  P3 
WRITE( 10,25)  T2 ,  P3 
WRITE( 10,27)  IP 
STOP 

20  FORMAT ( '  ENTER  THE  MAXIMUM  ALLOWABLE  TEMPERATURE,  ELSE  ' ,F7 .2, 'K' ) 

21  FORMAT ( '  THE  DIAMETER  OF  THE  LAYER,  ELSE  ',F7.2,'  cm') 

22  FORMAT ( '  THE  TIME  THE  BEAM  IS  ON  TARGET,  ELSE  ',F7.2,'  sec ' ) 

23  FORMAT ( '  THE  TOTAL  ABSORPTION,  ELSE  ',E7.1) 

24  FORMAT ( '  THE  LAYER  THICKNESS,  ELSE  ',F7.2,'  mm’) 

25  FORMAT ( /, '  THE  CENTRAL  TEMPERATURE  IS:  ',F7.2,'K', 

&  //,'  THE  MAXIMUM  POWER  DENSITY  IS  :'F9.4,'  MW/cm2 ' ) 

26  FORMAT(/, 2 OX, 'TEMPERATURE  PROFILE ’,//, 8X, < IR2 >X ,' RADI I  (cm.)'/ 

&  /,'  TIME  (  s  . )  '  ,  8  (  F7  .  2 , 2X) ,  /  ,  '  . ',<IRM>(2X,' . '/2X)) 

27  FORMAT ( /, '  IT  TOOK  ',14,'  STEPS  TO  FIND  THE  POWER') 

28  FORMAT( '  THE  SUM  LOOP  EXCEEDED  9999  TERMS') 

29  FORMAT)'  THE  TIME  INCREMENT,  ELSE  ',F7.4,'  sec') 

31  FORMAT ( '  RADIATION;  ',F9.2,'  EXCEEDED  TEMPERATURE;  ',F9.2) 

32  FORMAT ( '  THE  DIAMETER  OF  THE  BEAM,  ELSE  ',F7.3,'  cm') 

33  FORMAT ( '  THE  POWER  DENSITY  ',F9.2,'  KW,  GAVE  A  TEMPERATURE 
(.  '  ,F7 .2, 'K,  '  ,/,  '  WHICH  WAS  LESS  THAN  THE  MAX.  TEMPERATURE 
C.  '  ,  F7 . 2  ,  'K') 

91  FORMAT!'  THE  ABSORPTION  COEFFICIENT,  (ELSE  ',E7.1,'  per  cm)') 

92  FORMAT)'  THE  LAYER  THICKNESS,  IN  mm  (ELSE  ’‘IDEAL")') 

93  FORMAT)'  THE  PARTICLE  SI2E,  (ELSE  ',F5.2,'  urn)') 

94  FORMAT ('  THE  WAVELENGTH,  (ELSE  ' , F7 . 4 , '  urn)') 

95  FORMAT!'  THE  %-AGE  OF  SOLID  MATERIAL,  (ELSE  ',F3.0,'  %)') 

96  FORMAT( '  THE  INDEX  OF  REFRACTION,  (ELSE  *,F5.3,')') 

97  FORMAT)'  THE  MAXIMUM  ALLOWABLE  TRANSMISSION  *age,  IF  ANY.') 

98  FORMAT!’  FOR  Xi"  ',F7.3,'mm,  ',F8.5,'  l  IS  TRANSMITTED,  AND  ', 

(.  F8 . 5 , '  %  IS  ABSORBED') 

99  FORMAT!'  FOR  X-  ',F7.3,'mm,  ',F8.5,'  %  IS  TRANSMITTED,  AND  ', 

&  F8.5, '  %  IS  ABSORBED' ) 

120  FORMAT! 3F25. 16) 

121  FORMAT! 1X,F7.3,8(1X,F8.1)  ) 

END 
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o  n  o  *  n  n  r  n  n 


THIS  SUBROUTINE  CALCULATES  Jo( xB) /jo(B) **2 

FUNCTION  JO ( I ,  XX) 

IMPLICIT  REAL*8  (A-H,J,0-Z) 

REAL* 8  B( 10000) ,  J ( 99 9 9  ) ,  S(2) 

LOGICAL  OLD 
COMMON/  BETA  /  B 
COMMON/  LOGIC  /OLD 

ZRO*  0 .DO 
ONE*  1 . DO 
TWO*  2. DO 

PI*  3. 14 159 2653589793 DO 
P04*  PI  /  4. DO 

X=  XX  *  B(I) 

DO  30  IX*  1 , 2 

APPROXIMATION  FOR  LARGE  ARGUMENTS 

IF(  X  .LT.  6.D1  )  GO  TO  5 
S(IX)«  DSQRT(  TWO  /  (PI*X)  ) 

IF (  IX  .EQ.  1  )  S(IX)=  S(IX)  *  DCOS (  X  -  P04  ) 

GO  TO  30 

.  REVERSE  ITERATION  FOR  SMALL  ARGUMENTS 
C 

5  SUM*  ZRO 

IMAX=  INT( 9+1 . 83*X ) 

J ( IMAX+2 ) =  ZRO 
J( IMAX+1 ) =  l.D-30 
C 

DO  10  IM=IMAX,1,-1 

J  (  IM)  *  TWO*  IM*  J  (  IM+1 ) /X  -  J(IM*-2) 

IF ( IM  .EQ.  1)  GO  TO  20 
10  SUM*  SUM  +  (1  ♦  (-1)**( IM-1) )  *  J ( IM ) 

C 

20  SUM*  SUM  +  J( 1 ) 

S(IX)"  J ( 1 )  /  SUM 

X*  B(I) 

30  CONTINUE 
C 

J0=  S(l)  /  ( S( 2  )  *S ( 2  )  ) 

C  IF (  .NOT.  OLD  )  WRI TE ( 6 , 100 ) 

OLD-  .TRUE. 

C  WRITE (6,101)  I,  B ( I ) ,  S(l),  S(2),  JO 
RETURN 

100  FORMAT ( '  I ' , 8X, ’B(I) ',11X, 'Jo(rB) ' ,10X, 'Jo(B) ',12X, 'JO') 

101  FORMAT (lXfI4flX;4Fl6.10) 

END 


C 
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IDT=  MAX(  ITM/50  ,  1  ) 

DO  450  IT- 11, ITM 

IF(  IT  .GT.  ITM  )  GO  TO  500 

IF(  MOD( IT, IDT)  .NE.  0  )  GO  TO  450 

WRITE( 10, 121)  TIME ( IT ) ,  (T(IR,IT),  IR=1,IRM) 

450  CONTINUE 
C 

500  WRITE (6,25)  T2,  P3 
WRITE (10,25)  T2,  P3 
WRITE (10,27)  IP 
STOP 

20  FORMAT ( '  ENTER  THE  MAXIMUM  ALLOWABLE  TEMPERATURE,  ELSE  ',F7.2,'K') 

21  FORMAT (  '  THE  DIAMETER  OF  THE  LAYER,  ELSE  ',F7.2,'  cm') 

22  FORMAT ( '  THE  TIME  THE  BEAM  IS  ON  TARGET,  ELSE  ',F7.2,'  sec') 

23  FORMAT ( '  THE  TOTAL  ABSORPTION,  ELSE  ',E7.1) 

24  FORMAT ( '  THE  LAYER  THICKNESS,  ELSE  ',F7.2,'  mm') 

25  FORMAT ( / , '  THE  CENTRAL  TEMPERATURE  IS:  ',F7.2,'K', 

&  //,'  THE  MAXIMUM  POWER  DENSITY  IS  :'F9.4,'  MW/cm2 ' ) 

26  FORMAT(/,20X, 'TEMPERATURE  PROFILE ',//, 8X, <1 R2 >X ,' RADII  (cm.)', 

fi,  /,'  TIME  (  s  .  )  ' ,  8  (  F7 . 2 , 2X )  ,  /  ,  '  - ',<IRM>(2X,' - ',2X)) 

27  FORMAT ( / , *  IT  TOOK  ',14,'  STEPS  TO  FIND  THE  POWER') 

28  FORMAT { '  THE  SUM  LOOP  EXCEEDED  9999  TERMS') 

29  FORMAT ( '  THE  TIME  INCREMENT,  ELSE  ',F7.4,'  sec') 

31  FORMAT ( '  RADIATION:  ',F9.2,'  EXCEEDED  TEMPERATURE :  ',F9.2) 

32  FORMAT ( '  THE  DIAMETER  OF  THE  BEAM,  ELSE  ',F7.3,'  cm') 

33  FORMAT ( '  THE  POWER  DENSITY  ',F9.2,'  KW,  GAVE  A  TEMPERATURE 
&  ' , F7 . 2 , ' K,  '  ,  /  ,  '  WHICH  WAS  LESS  THAN  THE  MAX.  TEMPERATURE 
&  ' , F7 . 2 , 'K ' ) 

91  FORMAT ( '  THE  ABSORPTION  COEFFICIENT,  (ELSE  ',E7.1,'  per  cm)') 

92  FORMAT ( '  THE  LAYER  THICKNESS,  IN  mm  (ELSE  "IDEAL”)') 

93  FORMAT ( '  THE  PARTICLE  SIZE,  (ELSE  ',F5.2,'  urn) ' ) 

94  FORMAT ( '  THE  WAVELENGTH,  (ELSE  ',F7.4,'  urn)') 

95  FORMAT ( ’  THE  %-AGE  OF  SOLID  MATERIAL,  (ELSE  ',F3.0,'  %)') 

96  FORMAT ('  THE  INDEX  OF  REFRACTION,  (ELSE  ',F5.3,')') 

97  FORMAT ( '  THE  MAXIMUM  ALLOWABLE  TRANSMISSION  %age,  IF  ANY.') 

98  FORMAT ( '  FOR  Xi=  ',F7.3,'mm,  ',F8.5,'  %  IS  TRANSMITTED,  AND  ', 

fi,  F8.5, '  \  IS  ABSORBED' ) 

99  FORMAT ( '  FOR  X-  ',F7.3,'mm,  ',F8.5,'  %  IS  TRANSMITTED,  AND  ', 

&  F8.5, '  %  IS  ABSORBED' ) 

120  FORMAT ( 3F25-16) 

121  FORMAT ( 1X,F7.3,8(1X,F8.1)  ) 

END 

C 
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n  n  n  <r.  n  n  , ,  n  o 


THIS  SUBROUTINE  CALCULATES  Jo ( xB ) / Jo ( B ) *  *  2 

FUNCTION  J0( I ,  XX) 

IMPLICIT  REAL* 8  (A-H,J,0-Z) 

REAL*8  B( 10000) ,  J (  9 9 9 9  ) ,  S ( 2 ) 

LOGICAL  OLD 
COMMON/  BETA  /  B 
COMMON/  LOGIC  /  OLD 

ZRO=  O.'DO 
ONE“  1 . DO 
TWO=  2. DO 

PI-  3 . 14 159 26 5358979 3D0 
P04  =  PI  /  4. DO 

X*  XX  *  B ( I ) 

DO  30  IX-1 , 2 

APPROXIMATION  FOR  LARGE  ARGUMENTS 

IF(  X  .LT.  6 . DI  )  GO  TO  5 
S(IX)=  DSQRT(  TWO  /  (PI*X)  ) 

IF(  IX  .EQ.  1  )  S(IX)=  S (  IX  >  *  DCOS (  X  -  P04  ) 

GO  TO  30  , 

..  REVERSE  ITERATION  FOR  SMALL  ARGUMENTS 

C 

5  SUM=  ZRO 

IMAX=  INT(9+1.83*X) 

J ( IMAX+2 ) =  ZRO 
J ( IMAX+ 1 ) =  l.D-30 

0 

DO  10  IM=IMAX, 1,-1 

J(IM)=  TWO* IM*J(IM+1)/X  -  J ( IM+2 ) 

IF( IM  .EQ.  1}  GO  TO  20 
10  SUM=  SUM  *  (1  +  ( - 1 ) *  * ( IM- 1 ) )  *  J ( IM } 

20  SUM=  SUM  ♦  J(l) 

S(IX)=  J(l)  /  SUM 

X=  B ( I } 

30  CONTINUE 

JO-  S ( 1 )  /  (S { 2 ) *S ( 2  ) ) 

J  IF(  .NOT.  OLD  )  WRITE ( 6 , 100 ) 

OLD-  .TRUE. 

C  WRITE (6,101)  I,  B( I ) ,  S(l),  S ( 2 ) ,  JO 
RETURN 

100  FORMAT ( '  I' , 8X, 'B( I) ' , 11X, 'Jo(rB) ' , 10X, ' JO(B) ' , 12X, ' JO* ) 

101  FORMAT (1X,I4,1X,4F16.10) 

END 
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o  o 


THIS  SUBROUTINE  ENTERS  THE  FIRST  18  ZEROS  OF  THE  FIRST  ORDER  BESSEL 
FUNCTION 

C  FROM  Abramowitz  and  Stegun  "HANDBOOK  OF  MATHEMATICAL  FUNCTIONS"  AND 
APPROXIMATES  THE  REST 
C 


SUBROUTINE  BETAO 

REAL'S 

B(10000),  PI 

COMMON/ 

BETA  /  B 

PI*  3 . 

141592653589793D0 

B(l)  — 

3.8317059702 

B  (  2  )  = 

7.0155866698 

B  (  3 )  = 

10.1734681351 

B(  4  )  - 

13.3236919363 

B(  5 )  = 

16.4706300509 

B(  6 )  ® 

19.6158585105 

B  (  7  )  = 

22.7600843806 

B(8)  = 

25.9036720876 

B(9)  = 

29.0468285349 

B( 10 )  = 

32 . 1896799110 

B ( 11)= 

35.3323075501 

B( 12  )  = 

38.4747662348 

B( 13)= 

41.6170942128 

B ( 14  )  = 

44.7593189977 

B( 15)= 

47.9014608872 

B ( 16)= 

51.0435351836 

B  (  1 7  )  = 

54 .1855536411 

B( 18) = 

57.3275254379 

DO  1  1  = 

19,9999 

B(I)  = 

PI  *  (I  +  . 2 5D0 ) 

RETURN 

END 
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